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ABSTRACT 

We investigated the time-dependent radiative and dynamical properties of light su- 
personic jets launched into an external medium, using hydrodynamic simulations and 
numerical radiative transfer calculations. These involved various structural models for 
the ambient media, with density profiles appropriate for galactic and extragalactic 
systems. The radiative transfer formulation took full account of emission, absorption, 
re-emission, Faraday rotation and Faraday conversion explicitly. High time-resolution 
intensity maps were generated, frame- by- frame, to track the spatial hydrodynamical 
and radiative properties of the evolving jets. Intensity light curves were computed via 
integrating spatially over the emission maps. We apply the models to jets in active 
galactic nuclei (AGN). From the jet simulations and the time-dependent emission cal- 
culations we derived empirical relations for the emission intensity and size for jets at 
various evolutionary stages. The temporal properties of jet emission are not solely con- 
sequences of intrinsic variations in the hydrodynamics and thermal properties of the 
jet. They also depend on the interaction between the jet and the ambient medium. The 
interpretation of radio jet morphology therefore needs to take account of environmen- 
tal factors. Our calculations have also shown that the environmental interactions can 
affect specific emitting features, such as internal shocks and hotspots. Quantification of 
the temporal evolution and spatial distribution of these bright features, together with 
the derived relations between jet size and emission, would enable us to set constraints 
on the hydrodynamics of AGN and the structure of the ambient medium. 

Key words: radiative transfer — galaxies: active — galaxies: jets — quasars: 
general — radio continuum: galaxies 



1 INTRODUCTION 

Jets are present in astroph ysical systems from active galaxies 
(see iBegelman et al.l 1984ll to stellar systems, such as X-ray 
bin aries (Hiellming JohnstonI Il98ll : [Mirabel & Rodriguez? 
1999: Fender 2006) a nd symbiotic stars (Taylor et al. 1986; 
Sokoloski et^l200J ). They are supersonic. The jets from 



active galactic nuclei (AGN) are also relativistic, and they 
show a large variety of morphologies and properties. 

Although the jet plasmas have densities orders of mag- 
nitude smaller than the density of the ambient media, the 
large momentum flux carried in the plasmas enable the jets 
to propagate great distances, even up to Mpc scales. These 
jets impart momentum to the ambient gas. A termination 
shock (which appears as a hotspot) is often formed at the 
jet front, but in some cases, a series of lesser shocks along its 
length. Shocked jet plasmas fill low- density central cocoons, 
comprising eddies, backflows as well as transient shocks. 



Internal shocks may stand or propagate along th e jet it- 
self, c orresponding to intrinsic pinchin g; modes (e.g. Harded 
19791: ICohnlll983l: iNorman et al.lll982l: iHarded Il983l . Il987l : 



Falle et al.lll987l : lKoessl Muellerlll988l ) or stochastic dis- 
turbances from motions of the cocoon plasma. Studies have 
show that the termination shock may be stable or unsta- 
ble, depending on the system conditions. There are cases 
that the termination shock is suscept ible to vortex shed- 
ding, throbbing and disintegration (e.g.lNorman et al. 1982, 



Chakrabartilll988l : lKoessl Muellerll 19881 : iFalle 
i) 



1991 



Tf— 

O'Neill et al.ll2005l ). The disrupted hotspots and shocks in 
the backflow may be the cause of the bright rings or fila- 
ments seen in some radio galaxies. 

The energetic particles in jets emit strong synchrotron 
radio waves, and in some c ases, synchrotron X-rays (e.g. 
M87, [Marshall et aD l2002l : IPerlman WilsonI l2005l ). In 
FR2 radio galaxies, jets tend to end in termination shocks, 
where the jet plasma decelerates abruptly as it encounters a 
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dense external medium. The shocks compress the magnetic 
filed and convert the kinetic energy in the bulk flow into 
microscopic kinetic energy. In shock regions, charged 
particles are accelerated to relativist ic speeds, giving rise 
to synchrotron emission. The energetic particles in the 
jet also upscatter the low-energy photons (from the jet 
itself or from the ambient radiation field) to the X-ray 
energy ranges. Radio synch rotron emissions from j et s are 
linearly polarized (see e.g. iMantovani et al.l Il997l ). Some 
AGN also show substantial circular polari zation, in their 
core r adio emission and i n the j et knots (iHodge Ailed 
19771: IWeiler de Pate^ Il983l: iKomesaroff et all Il984l: 



Homan L Wardld Il999l: iRavner et al.l 12000171 Aller et al.l 
20031: iHoman Wardld l2004l: iHoman Listed l2006l : 



Vitrishchak et al.1 l2008l: ICenacchi et al.1 l2009l ). Theoretical 
studies (see McNamara et al.l l2009l ) have shown that 



Compton scattered X-rays from jets can also be strongly 
polarized. 

Jets are dynamical: their structures change over time. 
The jet evolution is manifested in morphological changes in 
their radio image and in their emission light curves. Stud- 
ies on the emission properties of jets tend to adopt sta- 
tionary structures. Simultaneous time-dependent radiative 
transfer calculations along with hydrodynamical simulations 
of jet evolution are rare as they are usually computation- 
ally demanding. However, only such time-dependent calcu- 
lations/simulations could allow us to generate appropriate 
mult i- wavelength lightcurve templates so to have a holistic 
pictures of the time-evolution of various jets (in particular, 
those in micro-quasars) and how they interact with the am- 
bient media. 

Here we present high time-resolution simultaneous ra- 
diative transfer calculations and hydrodynamic simulations 
of evolving jets in a variety of ambient media. The radiative 
transfer calculations are carried out frame by frame following 
the evolution of the jet density, velocity and thermal struc- 
tures in the hydrodynamic simulations. Time-sequences of 
emission images and light curves of jets are computed. Our 
study will shed light on the time-dependent hydrodynamics 
and radiative properties of radio jets. 

The paper is organised as follows. An introduction is 
given in §1; the methodology of the hydrodynamic simula- 
tions and radiative transfer calculations are laid out in §2. 
The emission properties of simulated jets, their interactions 
with ambient media, and the astrophysical implications are 
discussed in S3. We conclude in S4. 



2 CALCULATIONS 

2.1 hydrodynamic simulations 

In this study, we consider supersonic jets with low den- 
sity. Without losing generality, we consider axially symmet- 
ric jets. (Note that we have carried out full 3D hydrody- 
namic simulations for a few test cases without imposing 
the axi-symmetric condition. We have found that the re- 
sults are quantitatively similar to the corresponding axi- 
symmetric cases, but the computational time for each run 
does not let us conduct a comprehensive survey of various 
jet conditions and jet-interaction with the ambient media.) 
There are substantial works on the properties of such jets. 



fe.glNorman et al."l982': 'P e Youndll986l: ISmith et al.lll985l: 
Rosen et al.] [19 99: Carvalh ~ O'Deal l2002l : ISaxton et all 
2002ll?l: lKrau"s3l200a l2005|) • 

Our hydrodynamic simulations adopt the piecewise 
parabolic method (PPM) algorithm: an explicitly con- 
servative, shock-capturing;, grid-based numerical sc heme 
(jColella Woodward! Il984l: iBlondin LufkinI Il993ll We 



the P PMLR code ([Sutherland et al.ll2003l : ISaxton et al. lioOll . 
I2OO5I ) a descendant of the standard VH-1 University of Vir- 
ginia hydrocode[3 modified to improve numerical stability 
and parallelisation and to consider better tracer variables 
to distinguish fluids of differing origins. The computational 
grid consists of 600 x 300 cells (in longitudinal and ra- 
dial directions respectively). A reflecting condition applies 
on the cylindrical axis. Numerically, this border copies the 
nearby cells into corresponding ghost cells beyond the edge 
(with vector quantities reversed to ensure balanced restoring 
forces) . The right boundary and the outer cylindrical bound- 
ary are open to outflow. The jet nozzle spans 15 cells radially 
at the left boundary. Unless otherwise stated, a constant in- 
flow condition is applied. Outside the nozzle, the rest of the 
left boundary has either an "open" condition or a reflecting 
condition ("closed boundary"). For "open" boundary con- 
dition, we consider selective reflection to prevent inflow (if 
> at a border cell) . For edge outflow (v < ai the edge) , 
zero gradients (and zero force) is enforced, through copy- 
ing the single boundary cell into all ghost cells. No explicit 
treatment of radiative cooling is in the hydrodynamics for- 
mulation. We also omit self- gravity effects. These omissions 
would not cause significant effects and our results generally 
hold, provided that the free-fall timescales are longer than 
the crossing times associated with supersonic motions. 

In our formulation, the simulated jets are character- 
ized by two parameters: the jet Mach number M = v^/cs,} 
and the density contrast 77 = Pj/picm- Here Csj is the adia- 
batic sound speed in the jet plasma, and the subscript "icm" 
denotes the ambient medium. For adiabatic jets, variables 
such as mass, length and time can be rescaled, but the 
contrast variables 77 and A4) are fixed parameters. Gener- 
ally, the thermal velocity dispersion in the ambient medium 
may be used to set the velocity unit of the simulations. In 
this study, we consider an astrophysical realization appro- 
priate for giant radio galaxies. The cylindrical region in the 
simulations spans 150kpc x 75kpc. The jet radius is about 
3.75kpc. If we assume that the jet propagates in an intra- 
cluster medium (ICM) with a temperature of 10^ K and a 
density of 2 x 10~^^ g cm~^ (corresponding to a parti- 



cle density of nu 



10 cm ), then the velocity unit 



Uv = 365.5 km s~^ and the time unit Ut — 10.03 Myr. 

As a basis of comparison, we first carry out simulations 
of jets launched into uniform media. This simple setting 
avoids the unnecessary complications caused by the struc- 
tures in the ambient medium and enables us to identify the 
dynamics and emission properties intrinsic to the jets. We 
consider two jet models. In the first simulation (fastIo) 
we set = 50 and 77 = 10~^. The corresponding parti- 
cle column density through the jet is ^ 2.3 x 10^^ cm~^. 
In this case, the jets are expected to produce hotspots and 
other shocks with high brightness contrast relative to their 
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Table 1. Settings and parameters of the simulations. For those 
marked • the left boundary is closed/reflecting; for those marked 
O it is open. The jet density contrast and Mach number are 77 = 
Pj/Pism cind Ai. The duration of each simulation is tend cind the 
interval between frames is At (both in units of Ut ~ 10.03 Myr). 
Runs with nonuniform ambient media have index ^ and scale 
radius ro. (The grid size is 40 x 20 in the same units.) counts 
the frames. 
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surro undings (e.g. the radio galaxy Pictor A, ISaxton et alj 
I2OO2I ). In the second simulation (slow) we consider a lower 
Mach number [M — 3) for the jet. The density of the am- 
bient medium is the same. These jets tend to give lower 
luminous contrast for their induced f eatures (e.g. the radio 
galaxy Hercules A lSaxton et al.ir2002l V The upper and mid- 
dle rows of Figure [T]show density and pressure slices of SLOW 
and fastIo at late times. 

Next we run simulations for jets propagating into me- 
dia with various density structures. In the fast2 and fast2o 
simulations, we assume a decreasing medium density follo w- 
ing a flat-cored model (| Cavalier e Fusco-Femianolll976l ) 



1 



[l + (r/ro)2]3^/2 . 

We set the parameter ^ = 0.5 and the scale radius ro = 
15 kpc. Cases with a reflecting and an open left boundary 
are investigated. (Open cases are named with a suffix "o".) 
In the fast3 simulations, we consider a model where the 
density of the ambient medium increases with radius. We 
set ^ = — 1 and ro = 15 kpc. In the fast4 simulations, 
the density of the ambient medium varies sinusoidally, with 
p = po exp[— a sin(/cr)], where the wavenumber k = 27v/Xo. 
We set the amplitude a — log^ 10. This gives an effective 
wavelength of the variation Ao = 30 kpc. In the fastIr 
and fastIro simulations, we investigate the effects of the 
equation state. The jet plasma now has an adiabatic index 
7 = 4/3 (characteristics of relativist ic gas) instead of 5/3. 
Table [1] lists the principal parameters of the simulations and 
the temporal sequence of output frames. 

2.2 numerical radiative transfer 

In the radiative transfer calculations, we use a 3D Cartesian 
grid, consisting of 500 x 500 x 500 cells. The values of rel- 
evant variables on the radiative transfer grid are computed 
from the 2D data of axially symmetric jets obtained from 



the hydrodynamic simulations. For each Cartesian cell cen- 
tered at (x,^,z), we find the four neighbouring points of 
the cylindrical mesh that have integer coordinates. We use 
interpolation weights Wi oc l/{ri + 5) where Vi is the dis- 
tance neighbouring cylindrical vertex i and 5 = 10~^ is a 
softenning term in case the projected position lies near a 
cylindrical cell boundary or corner. 

The jets in this study are not magnetically dominated, 
and our simulations did not involve an explicit treatment of 
the magnetic field. However, the calculation of synchrotron 
radiation requires an explicit specification for the magnetic 
field. We therefore use a parametric prescription for the local 
magnetic field in the jet plasma. The field strength is taken 
as B = ySvrP/^, where P is the pressure and (3 (^ 1) is 
the strength parameter. The orientation of the field vector 
assumes an ordered component and a random component. 
The random component has no particular correlated length 
scales, i.e. with a white power density spectrum. 

The polarized radiative transfer is in the 4- Stokes for- 
mulation, with the Stokes parameters /, Q, [/ and V. The 
linear polarization component is given by L = y^W^ + Q^, 
and the circular polarization is given by V. The degrees 
of linear and circular polarization are Hl = L/I — 
yJTP^^^^~Q^ 1 1 and He — V/I respectively, where / is the 
total intensity. The Stokes parameters satisfy the time- 
dependent local polarized radiative transfer equations: 
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where , , and Sy are the corresponding emission coef- 
flcients of the Stokes parameters. The propagation operator 
is deflned as 
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where c is the speed of light, n is the normal vector of the 
ray propagation, I is the identity matrix, ds is the projected 
differentiation operator. The transfer matrix is given by 



(4) 



where k, is the absorption coefficient, q and v are coeffi- 
cients describing absorption effects; while / and h describe 
the Faraday propagation effects. The functional form of the 
emission, ab sorption, conver si on an d rotation coefficients 
are given m I J ones O'Delll (Il97/t ). Several studies have 
directly integrated the transfer equations (|2|)-(|4]) including 
circular polarisation, in parametri c models of shocks travers- 
ing cy lindrical or conical jets fe.g |jonesl[T988l : iHughes et aD 
Il989al lbl). We apply a similar treatment to hydrodynamically 
"live" simulated jets, but with a mathematical alteration to 
improve computational speed and robustness in inhomoge- 
neous structures of any opacity. 

To solve the transfer equation we use the following 
method. Consider the transformation: 

(d + K^I = J ^ \r~\D + K)r] R-^I = R-^J , (5) 

where / = [I Q U V]^ and J = [Sj Sq • Choose the 
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transformation matrix R such that it is translationally in- 
variant in space and time. The commutation ensures that the 
propagation operator D remains diagonahsed and is hence 
is invariant under the transformation. It follows that 



(d + k). 



(6) 



where / = R"^/, J = R"V, and K = R-^KR. When K is 
diagonalised, the components in the polarized transfer equa- 
tion are decoupled. Once the initial condition (the Stokes 
parameters of the incident ray) is specified, the equations 
can be integrated independently. This method is more com- 
putationally efficient and numerically stable than the direct 
integration of the coupled polarized radiative transfer equa- 
tion. 

Follow ing Wegd (l2003l), we use the algebraic package 
LAPACK (lAnderson et al.lll990l : fearker et al.ll200lh for di- 
agonalisation, which determines the eigen-matrix K and the 
transformation matrix R. Direct integration of Equation (|6]) 
yields /, of which an inverse transform under R gives the 
Stokes vector I. We consider a forward ray- tracing scheme, 
in which the transfer equation is solved sequentially along 
the ray through the cells. The collection of rays through the 
jet generates the time sequence of Stokes images. 



3 RESULTS AND DISCUSSIONS 

3.1 jet morphology and emission images 

3.1.1 standard cases 

Figure [2] shows the emission images of the fastIo and SLOW 
jets, calculated during advanced evolutionary stages. These 
two model jets are light and unobstructed. Bremsstrahlung 
maps are calculated by line-of-sight integration of the bolo- 
metric emissivity oc p^T^^^. In both simulations, most of the 
bremsstrahlung X-rays are emitted from the compressed am- 
bient medium within the bow-shock; the contribution of the 
jet itself to the total bremsstrahlung emission is insignifi- 
cant. The radio emission from the jets is synchrotron emis- 
sion, mostly from the hotspot, the jet's recollimation, pinch 
mode shocks, and the weaker shocks in the turbulent co- 
coon. In producing 6 cm radio synchrotron emission maps 
of these jets, we adopt a relativist ic electron energy spectral 
index p = 2.2. By default we have assumed /3 = 10 (i.e. 
hydrostatic pressure dominates magnetic pressure) . In some 
extra comparisons, we also calculate the maximal field case 
{/3 = 1) or a weaker field case (/3 = 100). 

F or jet parameters equivalent to fastIo, I Saxton et al.l 
(I2OO2I ) noted axisymmetric shocks spawned off the throb- 
bing hotspot, and transverse filaments emerged, resembling 
structures seen in the western lobe of the radio galaxy Pic- 
tor A. Here we reproduce comparable models as a bench- 
mark. In addition, we have carried out explicit numerical 
polarized radiative transfer calculations frame by frame to 
show the emission and structural evolution of the jets. Fig- 
ure [3] shows a sample sequence of image frames of the radio 
intensities of the evolving fastIo jet. The morphology is 
relatively simple at the early times, with the emission domi- 
nated by the hotspot. As the jet evolves and propagates, in- 
ternal recollimation shocks begin to appear, and filamentary 
shocks are thrown off the jet's advancing head. At certain 
epochs, the filamentary shocks appear ahead of the hotspot. 



but at other times they are coincident or behind (resem- 
bling structures in Pictor A. iRoeser Meisenheime 3 ll987l : 
IPerlev et al.lll997l : IWilson et al.ll200ll : iTingav et al.1 120081 ). 

The hotspot and recollimation shocks appear sharp in con- 
trast to the cocoon. Reassuringly, the synchrotron morpholo- 
gies generated by full radiative t ransfer resemble the approx- 
imate imaging bv I Saxton et al ] §002). 

Radio intensity maps of the SLOW jets show more inter- 
nal shocks for a given length of jet. Globally the images are 
less contrasty than for the fast jets. There is a conspicuous, 
fuzzy bulge is located near the jet base. It is turbulent jet- 
derived plasma in the backflow around the jet. When the 
boundary is open, allowing outflow off the grid (e.g. com- 
pare the fast2o jets to the fast2 jets), the fuzzy feature 
diminishes. We suspect that it might diminish further in the 
3D simulations and when the dimming of aged plasma were 
taken into account. The hotspot 's local behaviours are ef- 
fectively independent of the left boundary condition when 
the jet has progressed far enough, in the late evolutionary 
epochs of each simulation. 

Our radiative transfer calculations show that the jet 
emission structure depends on the assumed model of 
the magnetic-field configurations. If we assign the local 
magnetic- field orientations parallel to the velocity vectors, 
the intensity distribution in the jet image will be quite 
smooth (after the subtraction of bright artefacts that occa- 
sionally appear in the cocoon) . If a random field orientation 
is used (with directions independent between adjacent cells) 
then the jet image will appear speckly. In ^3.1.3l we consider 
the effect of some alternative configurations of the ordered 
magnetic field. 

In some simulations, large eddies with high concentra- 
tions of jet plasma are spun into large rings (by axisym- 
metry). Projections of these structures gives rise to patches 
with intense emission. The similarity of the light-curves in 
the "random" and "ordered" cases of Figures [714151 demon- 
strates that these twinkles cause only local fluctuations. 
Their effects on spatially integrated intensity and Stokes pa- 
rameters are usually negligible. 

For synchrotron emission from relativistic electrons 
with a power-law energy distribution of an index p gyrat- 
ing in an orderly oriented magnetic field, the intensity spec- 
tral index a and the linear polarization IIl are given by 
a = (p — l)/2 and IIl = (p + l)/(p + |) respectively. As- 
suming p = 2.2 in the calculations gives an index a = 0.6 of 
the synchrotron spectra and a maximum linear polarization 
IIl ^ 0.7. The peak local values of IIl in the simulated jets 
generally approach the theoretical maximum value of 0.7 
regardless of the assigned magnetic- field configurations, as 
expected. Sy nchrotron emission has little intrinsic circular 
polarization (iLegg k Westfoldlll968l ). However in a plasma, 
some circular polarization can be produced by Faraday con- 
version of the linear polarization components of radiation. 
(We discuss an aspect of polarization evolution briefly in 
^3.41 and in more detail in a forthcoming paper.) 

Jet pinch shocks show a simple luminous knot where 
the biconical shock occurs. Hot-spots are composed of more 
complex and transient shocks, resulting in variable, closely 
spaced luminous patches. Resembling FRl radio galaxies, 
the SLOW jets do not have a clear hotspot, and so the co- 
coon dominates the radio emission. For random fields, the 
jet is inconspicuous. There is little cocoon obscuration in the 
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Figure 1. Sections of jet through the cyhndrical structure at final frames of three simulations. The grey values and scale bar are decimal 
logarithmic. The jet axis runs along the bottom of each frame. The radial direction runs upwards. The left panels show the density; the 
right panels show the pressure (a proxy for B^/Stt). The SLOW jet is on the upper, and the fastIo jet is on the middle row. The jet 
plasma is a polytropic gas with 7 = 6/8. The fastIro jet is on the bottom row. The simulation has the parameters and set ups same as 
those of the fastIo jet, except that the jet plasma is relativistic with 7 = 4/3. 



fastIo jets, and the jet shocks and the hotspot are distinct, 
like the morphology of FR2 sources. 



3.1.2 alternative environments 

We now comment on the characteristics of several variants 
of simulated fast jets: 

fastI: The jets propagate in a uniform background den- 
sity, as was the case for the fast lo jets. However a closed left 
boundary is now imposed. For FASTl jets, the jet itself and 
the shocks are generally distinguishable, but cocoon emis- 
sion is more apparent. Many pinch shocks appear in the jet 
at the late stages of the evolution. The hotspot surges irreg- 
ularly back and forth, typically traversing about a third the 
jet's length. 

fastIro: The background is uniform and the left 
boundary is open to outflow, but the adiabatic index 7 = 



4/3. Compared to fastIo, the cocoon is narrower and the 
jet takes longer to cross the grid. 

fastIr: This is equivalent to fastIro except that the 
left boundary is closed. Closure causes a wider cocoon to 
accumulate (as in fastI), despite the choice of 7 = 4/3. 

fast2: The jets propagate in a radially declining back- 
ground density. The jet emission remains distinctly visible 
through a cloudy cocoon (e.g. Figure |4]). There are fewer 
pinch shocks than the case with a uniform background den- 
sity, and the shocks are farther apart. The first two pinch 
shocks tend to outshine the hotspot. The hotspot throws off 
rings or filaments, as in the case with a uniform background 
density. 

FASt2o: An open boundary reduces the emission inten- 
sity of the cocoon. Apart from the reduction of the emis- 
sion of the cocoon, the emission structure of the fast2o jets 
resembles that of the fast2 jets. The jet and hotspot are 
clearer than those of the fast2 jets. 
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Figure 2. Emission images of the fastIo jet (top row) and the SLOW jet (bottom) at their late evolutionary stage. The image sizes are 
150 kpc X 150 kpc. In both cases, /3 = 10 and the field is quasi-poloidal (B \\ v). Columns from left to right show: thermal bremsstrahlung 
X-rays, and synchrotron radio intensity (/). The emissions are linearly scaled. The synchrotron intensity is evaluated at 6 cm. The jet 
morphology in the non-thermal synchrotron X-ray images (unshown) is similar to that in the synchrotron radio emission image. 



fast3: a rising background density confines the expan- 
sion of the turbulent cocoon. The jet's advance and cocoon's 
growth decelerate (the front of the bow shock Xs oc t^'^^). We 
stop the simulation before the jet traverses the grid. With 
ordered fields and /3 = 10, the late-time images only show a 
translucent cocoon with eddy induced substructures. (The 
case with /3 = 1 is unassessable, because the twinkle arte- 
facts dominate.) The jet is obscured and there is no apparent 
hotspot (see Figure [5]). 



FAST4: The background medium has a ripply density 
variation. Naturally, the intensity of the jet emission varies 
as the jet propagates into these shells. In spite of the density 
variation, the jet morphologies are similar to those of the 
fastI jets. The emission of the cocoon dominates in the 
intensity image, as in all cases with a closed left boundary. 
There is some subtle pinching of the cocoon outline where 
the jet reaches breakthrough points. At later evolutionary 
epochs, the jet head surges stochastically as far as half the 
jet's length. 



3.1.3 magnetic field configuration 

In our standard calculations, we have assumed that the mag- 
netic field is parallel to the flow, i.e. B \\ v. The field strength 
is parametrized, and its energy density scales with plasma 
thermal pressure. This corresponds to a quasi-poloidal field 
in the jet, but with "live" deviations at transient features 
such as shocks and vortices. As radio synchrotron radiation 
is determined by the magnetic field, we therefore investigate 
the jet emission for other magnetic field configurations. As 
a test of generality, we perform additional calculations with 
two alternative field configurations. (Figure [6|). 

In the sequence of the TURNl simulations, we rerun cal- 
culations for the fastIo simulations, with magnetic field 
directions perpendicular to local velocities. The cylindrical 
components are B = {Bz, Br, B^) \\ {vr, —Vz,0). Fields in 
the jet tend converge towards the axis (a quasi-radial con- 
figuration) which may be less physical than the usual B 
recipe. The emission morphology is essentially the same as 
that of the fastIo jets (top panels. Figure [6|), but with a 
dark midline due to symmetry in B. 

In a futher sequence of the curlI simulations, we set 
B II V X V. This prescription guarantees V ■ B — 0, and 
yields quasi-toroidal fields around the jet. In the intensity 
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Figure 3. A time sequence of synchrotron radio intensity images of the fastIo jet at various evolutionary stages (frames 54, 108, 161, 
215, 269, 323, 376 and 430, corresponding to times t ^ 0.217, 0.433, 0.646, 0.863, 1.08, 1.30, 1.51, 1.73 Myr after the initial launch of the 
jet). In this calculation, (3 = 10 and the field is ordered (quasi-poloidally B \\ v). 




Figure 4. The last frame of the FAST2 jet simulation. In the 
simulation, (3 = 10, and the local magnetic field orientation is 
aligned with the fiow velocity. Note that the hotspot is developing 
into a funnel morphology. 



Figure 5. The 6 cm synchrotron radio intensity (/) map of the 
final frame of the FAST3 jet simulation. In this simulation, the jet 
encounters a rising density profile. The magnetic field is quasi- 
poloidal {B \\ v) and (3 = 10. 



images, the brightness is concentrated closer to the axis, 
but knots and the hotspot still appear in the usual places 
(bottom panels, Figure [6]). 

The general appearance of the jet and local features do 
not vary significantly under drastic changes of the B config- 
uration. The 6 cm radio emission morphologies are limited 
by the underlying plasma density and pressure distributions. 
Within the luminous features, there are cosmetic variations 
caused by the minor polarisation effects and magnetic con- 
figuration. Mindful of these small differences, we proceed 
under the fiducial, quasi-poloidal assumption, and calculate 
radio maps and light-curves of fast jets penetrating various 
non-uniform media. 



3.2 Effects on jet evolution and emission 
variability 

Figure [3 shows the evolution of the integrated intensity 
/, of the SLOW jet. Figures [8lll5l are equivalent light-curves 
for the fast jets. In simulations where the hotspot is unam- 
biguous and unobscured by the cocoon, we also present a 
local light-curve of the hotspot (/hs(t)). Here we define the 
instantaneous "hotspot" as the intensity maximum that is 
farthest right along the jet axis, and we integrate the emis- 
sion within a projected radius of 20 pixels. 
3.2.1 left boundary conditions 

We have considered two types of left boundary con- 
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Figure 6. Synchrotron radio intensity (/) maps (at 6 cm) of the 
fast jet as in Figure[2] (fastIo), but with a different prescription 
for the magnetic field geometry. Top panel shows turnI, with a 
quasi-radial field (B ± v). Bottom panel shows the result of a 
quasi-toroidal recipe (CURLI, where B \\ V x v). 



ditions in the simulations. The open boundary condition 
mimics situations where the jets have drilled far from the 
central engine. This is appropriate for modeling the jets in 
giant radio galaxies. The closed boundary takes account of 
the effects of backflows and material accumulation. This is 
appropriate for modeling compact jets that are yet to break 
away from the core region, where backflows from both jet 
and counter-jet can interact with each other. Comparison of 
runs that differ at the left boundary (i.e. fastIo vs fastI 
and fast2o vs fast2) shows that the (spatially integrated) 
temporal fluctuations are more rapid and have larger am- 
plitudes for an open left boundary. This may be attributed 
to fact that the exit of plasma through the boundary leaves 
the unsteady structures (such as the jet shocks and hotspot) 
brighter than the cocoon. Also, with an open boundary con- 
dition, the cocoon contributes less to the total emission. 

The cloudy cocoons of the 7 = 5/3 and 7 = 4/3 cases 
with closed boundaries have similar light curves (see Fig- 
ures [9] and [11]). The time taken to traverse the grid is not sig- 
nificantly affected by the alternative 7. The intensity reaches 
marginally higher levels for 7 = 5/3. When the left boundary 
is open (Figures [8l and [TQ|) the effect of 7 is more conspicu- 



ous: 7 = 4/3 gives greater variability in the total intensity 
/(t), as well as in bright features like the hotspot Ihs{t). 

3.2.2 magnetic field structure 

The differences in the intensity light curves between a 
random or an ordered field directions are more obvious for 
jets with a closed boundary and conspicuous cocoon (e.g. 
comparing (5—10 columns within Figures [9l 1131 [T^ and 
[T5|) . Integrated intensity and polarization profiles tend to 
be smoother in the random case, due to summation over 
incoherent zones of the image. Ordered fields following the 
velocity structure produce larger coherent emssion patches 
that evolve on the same temporal and spatial scales as the 
principal shocks and eddies. The hydrodynamic vacillation 
and surging of these structures is responsible for the high 
frequency quasi-periodicity apparent in the I{t) curves. 

3.2.3 approximate homology of intensities 

A comparison of intensity light-curves calculated from the 
same simulation but with different /3 shows that they are 
similar to within a scale factor of ~ 0.16 for each factor 
of 10 increase in (5. For cases with a closed boundary and 
a translucent cocoon (slow, fastI, fast2, fast3, fast4, 
fastIr) the homology holds at nearly all times. For the 
cases with an open left boundary and prominent jet shocks, 
the homology holds at late times (say when the front of the 
bow shock is beyond Xs > 30 kpc) but the initial behaviours 
differ in detail. The scaling also breaks down temporarily 
around particularly bright and opaque features during the 
undulations of fast4b and flashes of the jet in fastIro. 
The extreme choice of /3 = 1 (strong field) makes more sub- 
structures opaque, and the scaling fails more often for higher 
13. 

3.3 Jet variability and astrophysical implications 

3.3.1 coruscating substructures 

Evolution of local substructures along the midline of the im- 
age is shown by intensity slices in Figure [16] The hotspot has 
rapid structural and morphological variabilities: it throbs; it 
casts off vortex rings; it breaks and reforms. Identification of 
the hotspot is not always unambiguous. In one case (fast3), 
the hotspot and jet internal shocks are hidden within an 
opaque fireball. In general, the jet shocks have comparable 
intensities for cases with an open left boundary. However, 
when when left boundary is closed (fastI, fast2, fast4) 
the most persistently bright feature is the first jet shock. It 
outshines the jet hotspot, a nd would appear as a "pseudo- 
core" as in jet observations (|jonesl[l988l : IP'Arcangelo et al.l 
I2OO7I ) . At later evolutionary stages, a second pinching shock 
appears. It could be as bright as the hotspot, and it tends 
to migrate forwards (whereas the first shock hovers about a 
certain place). 

Subtler diagonal streaks through the diagram represent 
more ephemeral shocks traversing the jet. The formation of 
bright emission knots in jets is sensitive to the environment, 
and the distribution of bright knots in the jets during any 
particular epoch could be a diagnostic of the density struc- 
ture of the external medium. Our simulations have shown 
that the knots occur at nearly equal intervals for a uniform 
background (e.g. fastI, fastIo). The knots can also appear 
evenly spaced, for jets propagating through a medium with 
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density that undulates about a mean (fast4). For a medium 
with a radially declining density (fast2o), the knots are 
closer together in the outskirts, beyond a long gap in the 
denser inner region. 

3.3.2 intensity and jet advance 

As shown by each wedge's outline in Figure [161 the radio- 
emitting region advances as time elapses, and the front of the 
bow shock (denoted as Xs, not illustrated) evolves similarly. 
Their evolutions are dependent on the external density pro- 
file. In a uniform medium, the jet advances and the cocoon 
expands with a gentle deceleration. The deceleration is more 
severe for media with a radially rising density (fast3). How- 
ever, acceleration would occur if the density of the ambient 
medium decreases radially (fast2, fast2o). 

In nature, the age of a jet activity episode is not di- 
rectly observable, and neither are the density profiles of the 
medium in the vicinity (say, within a few tens kpc) of an 
AGN core. In contrast, the size of the system and the ra- 
dio luminosities can be estimated from observational data. 
Thus, deriving a relation between the system size and ra- 
dio luminosity from our simulations of time-dependent jets 
in various environments will provide a tool to constrain the 
age of jet activity episodes, and infer the ambient structure. 

For models with an open left boundary, the light-curves 
are dominated by stochastic fiashes from jet shocks. Once 
such a jet has grown long enough to have at least one pinch 
shock, there isn't any tight relation between size and lumi- 
nosity. (For FASTlo, fast2o and fastIro the light- curve 
at late times brightens and dims significantly on diverse 
timescales.) However for the models with a closed boundary 
(forming a foggy cocoon/lobe) relatively simple power- law 
relations occur. 

The integrated intensity light curves I{t) of the fastI 
and fast2 jets are quite similar except that the latter 
drops at later times, when the jet encounters a signifi- 
cant density drop. During the initial rising stage of fastI 
(6.9 kpc < Xs < 15 kpc) the total intensity grows like the 
plasma volume / oc oc t^"^^. The instantaneous peak 

intensity varies rapidly within a broad envelope that tends 
to run like /max oc x^'^^ oc t~°"^^. The growth curve steep- 
ens over time: Xs oc t°"^^ early on (9 kpc < Xg < 15 kpc). 
At later times (19.5 kpc < Xs < 120 kpc) the growth is 
Xs oc t°"^^, while the luminosity / oc x^'^^ oc and the 
peak /max oc x~^-^^ oc t~^-^'^. 

For fast2, the Xs = Xs(t) curve steepens as the jet ad- 
vances, but the index is > 1 for most of the simulation. When 
Xs > 30 kpc the mean log-slope is 1.25; when Xs > 75 kpc 
this rises to 1.53. During the initial rise (xs < 12 kpc) 
the luminosity evolves as / oc xf ^^ oc t^"^^ and the peak 
/max oc Xg-^^ OC t^'^^ . After Xs > 12 kpc, the total inten- 
sity curve is nearly fiat; the instantaneous peak values drop 
roughly like /max oc x~^-^^. 

The fast3 jet infiates a globular cocoon/fireball with 
plasma accumulating at a constant rate, but a decelerating 
radial expansion. At later stages (xs > 27 kpc), the growth 
curve is Xs oc t°"^^. As the cocoon is opaque for /3 = 1, 10 or 
100, the luminosity rises with the surface area, / oc x^ *^^ oc 
^0.92^ (Figure [T4|) . The (noisy) envelope around the peak 
intensities rises slower, /max oc x'^'^^ oc t^'^^ . 

The FAST4 jets (Figure [T5|) . propagating in a ripply 



background, have light curves similar to fastI, but with un- 
dulations, as the jet is retarded by overdense shells, and then 
surges faster through underdense shells. 

The light curves of the SLOW jets resemble the fastI 
and fast2 jets initially and fast4 jets in the later stage. 
Despite the SLOW jet's low Mach number, it advances with 
a power- law behaviour resembling that of fastI: Xs oc t°"^^ 
at late stages. The intensity laws are steeper than for fastI, 
/ oc Xg'^^ and /max oc x^*^"^^- During the early rise, Xs oc 
/ oc a;2-26 oc f-^' and 1^,^ oc 

3.3.3 hotspot variability 

For models with an open boundary, the cocoon is mostly 
optically thin and the hotspot and other bright knots are 
relatively bare, so it is feasible to characterise their variabil- 
ities. In the top row of Figure [TTl we plot the rms variability 
of the hotspot intensity, calculated from the temporal struc- 
ture function of Ihs{t). 

&(/hs;r) = — ^ T'"^ [/hs(t + T)-/hs(t)]' dt. (7) 

tend T JQ 

The bottom row shows corresponding structure functions 
calculated from the total integrated intensity, /(t). 

The basic jet in a uniform background (fastIo) has a 
structure function that rises towards longer timescales. Per- 
haps this refiects the trend for later fiashes of the hotspots 
to be dimmer than earlier fiashes (row two of Figure [8]). 
The structure function is nonzero down to the smallest r 
scales; implying that the time steps of our hydrodynamic 
snapshots are too coarse to resolve the vacillations of the 
brightest shocks. Changing the field geometry from random 
to ordered (e.g. with /3 = 10 fixed) has no significant effect 
on the structure functions. Greater /3 (weaker fields) raises 
the size of the Ihs fiuct nations relative to the mean value, 
but the structure function retains the same shape in r. 

With a declining density profile (fast2o) the structure 
function of Ihs reveals peaks around timescales r ^ 0.05Ut 
and O.lGUt. The 0.05Ut bump may refiect the long in- 
terval between prolonged bright fiares (seen in the initial 
episode and two later events in Figure fT2)) . The peaks be- 
come sharper and more distinct for greater /3: given the same 
fiow evolution, a less magnetised hotspot fiashes relatively 
more intensely. The structure function of total intensity /(t) 
increases at longer timescales in the domain r < O.lOUt, per- 
haps describing with the brightenning trend, or the jump in 
/ around the time (~ O.OSUt) when the hotspot first sepa- 
rates from the pseudocore shock. The separation of subse- 
quent diamond shocks at ^ 0.02Ut intervals do not produce 
a clear feature in the structure function of /. 

The jet with a fiat background and relativist ic equation 
of state (fastIro) has a hotspot that fiuctuates with similar 
power at all time scales. The structure function of Ihs{t) is 
similarly fiat for /3 = 1, 10 and 100. There is a slight hump 
at nearly the duration of the simulation; this may be due to 
the early, dim and steady phase at t < 0.9Ut. The general 
fiatness of the temporal structure functions suggests that 
the stochastic variation of the hotspot is like white noise (at 
least for the "hotspot" as we have presently defined it). 

A radio galaxy of FR2 type has two opposite hotspots. 
Even if both the jet and the counter-jet are identical and un- 
varying at the nucleus, their hotspots vary independently: 
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differently at any given epoch, but stochastically in essen- 
tially the same manner. We assume that the envelope or 
distribution of this innate variability is well characterised 
by the recorded variability during sufficiently late stages of 
the simulations. Frames occuring after the formation of the 
pseudocore shock should be an adequate selection. 

The upper panels of Figure [18] show what is effectively 
the goemetric mean of the hotspot ratios at a time separa- 
tion T, as calculated from the structure function of the log- 
intensities, exp y^5'2(ln(/hs); t). Typical values are ^ 3 for 
fast2o and fastIro, but exceed ten at medium timescales 
for fastIo. In the lower panels of Figure [18] we draw ran- 
dom values of Ihs from times when Xs ^ 60 kpc, and plot 
the cumulative distribution function of the hotspot lumi- 
nosity ratio. Our radiative transfer calculations imply that 
ratios of factors of a few can occur frequently by chance 
(even though the paired jets are equal and steady). Ratios 
in excess of ten are not rare for fastIo. For each simulation, 
the 13 = 1 calculations give the least variable hotspot ratios. 
Perhaps this is because stronger magnetism raises the opac- 
ity of the brightest features, so that the foreground emitting 
surfaces effectively hide some of the variability farther from 
the virtual camera. The 7 = 4/3 simulation fastIro shows 
this effect most strongly. The model with a radially declin- 
ing density profile for the ambient medium (fast2o) shows 
less variable hotspot ratios than the standard uniform back- 
ground (fastIo). 

The general lesson to be drawn from these plots is that 
the apparent differences between opposite hotspots (during 
a single observational epoch) may be largely stochastic and 
transitory (even while the jet nozzles remain steady) and 
that caution is warranted when using these comparisons 
alone to constrain jet orientation, power or the behaviour of 
the nucleus. Doppler and orientation effects complicate mat- 
ters further, calling for more arduous parameter surveys in 
future. To make precise observational inferrences, it may be 
best to involve complementary evidence from other bands, 
or to take a statistical view of large ensembles of AGN. 

3.4 stages of jet activity 

Adiabatic hydrodynamic simulations are rescalable in the 
three physical units: length, time and density. The imple- 
mentation of radiative transfer calculations, however, lim- 
its this freedom, because of the constraint that homolo- 
gous models have the same optical depths (at correspond- 
ing scaled radiation frequencies) after rescaling the physical 
structure. Nonetheless, some freedom to rescale the models 
remains. Our models are therefore qualitatively relevant to 
analogous systems at smaller sizes, such as microquasar jets 
(although till this section we have interpreted our simula- 
tions with physical scales typical of large radio galaxies). 
For instance, a parsec-scale dwarf version of the confined 
fireball fast3 is conceivable. 

The earliest stage of jet expansion involves a compact 
ball of jet plasma, with at most a single vortex cell (a much 
simpler backflow than at late times). A prominent termi- 
nation shock shines brightly, but the jet is too short to 
fit any series of internal shocks. This is probably what oc- 
curs during flares by AGN and analogous microquasars (e.g. 
Fender et al. |200Q| . l2QQ2l : iMacailart et al I I2OO2I : lAller et al.1 
20031 ). As Figure 1161 shows, for the fast jets with an open 



boundary the terminal shock dominates the luminosity dur- 
ing this stage (t < 0.02Ut for fastIo; t < OmUt for 
FASt2o). In corresponding models with a closed boundary 
(fastI; FASt2) the fireball is more uniformly luminous due 
to the accumulated mass of plasma surrounding the sides of 
the jet. 

The closed-boundary simulations represent jets still 
propagating near the nucleus. The accumulated plasma co- 
coon contributes a fuzz of radio emission on lateral scales 
a considerable fraction of the jet's length. After the ini- 
tial fireball stage, the locally brightest feature is usually the 
pseudo-core (first jet internal shock) but the extended fuzz 
contributes most of the integrated emission. In fast4 the 
ripply density background modulates the cocoon expansion 
dramatically (and causes a long-period throbbing in the I{t) 
light curve), yet the position and variability of the pseudo- 
core are indistinguishable from the simpler fastI case. In 
integrated intensity observations of closed systems, the fuzzy 
cocoon may wash out the intrinsic variability of the jet (as 
seen in the slowly varying /(t) light curves). Spectral ag- 
ing and dimming of the cocoon plasma (not included in our 
present simulations) would leave the jet relatively prominent 
at later times. 

Simulations with an open left boundary represent jets 
that have penetrated far from the nucleus. Jet internal 
shocks and the hotspot then dominate the intensity maps 
at all times after the initial fireball. A comparison of fastI 
with fastIo or fast2 with fast2o shows a greater sep- 
aration between jet internal shocks in cases with an open 
boundary. When the background density profile decreases 
radially (fast2o) the pseudo-core is farther from the nucleus 
but the subsequent jet-shocks are closer together (compared 
to FASTlo). 

What about a jet injected into a pre-existing cavity, 
such as the ghost cocoon from a previous episode of nuclear 
activity? If the jet meets a smoothly rising density profile of 
the interstellar medium (like the case of fast3) a very differ- 
ent morphology emerges: a prolonged, highly opaque, frus- 
trated fireball. This spheroid is edge-dim, unlike an early- 
stage fireball with bright termination shock. The frustrated 
fireball contains a more complex system of eddies, with dark 
creases visible at eddy interfaces nearer the surface. If these 
flows carry locally ordered magnetic fields, then these affect 
the radiative transfer differently to the simple fields in an 
"early fireball" . The "frustrated old fireball" and "early fire- 
ball" may have distinct polarization signatures. (We study 
these in detail in a forthcoming paper.) 

As shown in Figure [TBI although the internal jet shocks 
may hover about preferred points, they eventually drift up- 
stream or downstream at apparently high velocities. Our in- 
tegrated light profiles and hotspot light-curves for the open- 
boundary fast jets (upper two rows of Figures [8] [10] and 
I12p show that internal shocks and hotspots can brighten or 
fade by factors of a few on timescales of millennia or less. In 
smaller analogues of AGN, with shorter timescales, it is con- 
ceivable that an internal shock that fades, moves upstream 
and rebrightens might be mistaken for a new ejection there. 

We noted that a pulse of circular polarization (He) can 
occur in the very initial stage of the simulations (bottom row 
of Figure [HJ, when the jet has a relatively coherent struc- 
ture. The occurrence of a circular pulse is robust and seems 
insensitive to the density profile of the external medium. It 
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Figure 7. Light curves of the simulated SLOW jet (uniform background; closed boundary). The left column show the case with a random 
field orientation and /3 = 10. The middle and right columns show cases with ordered local fields aligned with the matter fiow velocities, 
with (3 = 10, 1 respectively. The integrated 6 cm synchrotron radio intensity / light curves are on the top panels. Note that the fuzzy 
cocoon dominates the emission, masking variability on short timescales. 




Figure 8. Intensity and polarization light-curves from radiative transfer calculations based on the fastIo jet simulations, where a 
fast jet propagates in a uniform background medium. An open left boundary is adopted. The jet emission is dominated by distinct and 
coherent features, such as the hotspot and jet shocks. Dots mark exact data; grey curves show a 5-point boxcar average. The spatially 
integrated intensity light curve (first row) is variable on many time scales. Hotspot intensity (/hs? within a 20-cell radius) light curves 
are in the second row. The hotspot is identified as the rightmost local maximum of intensity projected along the jet axis. The /hs light 
curves fiuctuate considerably more than those of the SLOW jet and the fastI jet with closed boundary, as the hotspot of the fastIo jet is 
unobscured by the cocoon. The third and fourth rows show the integrated fractional linear and circular polarisation curves respectively. 



occurs whether the magnetic configuration is set to be quasi- 
poloidal (dragged parallel to the jet flow) or quasi-toroidal 
(circulating around the flow), but not if the fields are ran- 
domly directed. Later, the circular polarisation is diluted 
rapidly as the jet evolves more numerous, complex, mutu- 



ally incoherent emitting substructures. Thus, the circular 
polarizati on observed in A GN cores (e.g. iHoman Wardl3 
1999; Ra vner et al.l l2000l ) probably indicates certain de- 
grees of coherence in the emission region. Parametric mod- 
els with turbulent or ordered, poloidal, toroidal or (in- 
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Figure 9. Same as Figure [Jl but for the fastI jet simulations. In the simulations, a fast jet propagates in uniform background medium, 
and a closed boundary condition is adopted. Compared to fastIo, there is more emission from the cocoon. The lightcurves thus track 
the expansion and emission fluctuations of the cocoon rather than the variability in the jet shocks. 
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Figure 10. Same as Figure [71 but for the fastIro jet simulations. In the simulations, a fast jet propagate in a uniform background 
medium, and an open boundary is adopted. The jet plasma is relativistic, with 7 = 4/3. 



termediately) helical magnetic configurations have previ- 
ousl y been applied to explain polarised emissions of AGN 



fe.e. iBeckert & Falckell2002l:| 


luszkowski & Beeelman 2002"; 


EnBlinI l2003l: iGabuzda et alJ 


20081: IVitrishchak et al 2008i: 


Homan et alJ bood. and references therein). The essential 



ingredient in these models and in our time-dependent mod- 
els is that magnetic fields vary along sightlines through the 
jet, so that linearly polarised radiation from distant loca- 
tions undergoes Faraday conversion to circular modes in 
foreground plasma. As expected for this mechanism, we see 
that each spike in |nc| coincides with a dip in linear po- 
larisation IHlI (e.g. third row. Figure [8|). Circular polariza- 
tion at a few percentage level was also seen in the initial 
stage of radio outburs t in the micro quasar GRO J 1655—40 
(jMacQuart et al.ll2002l ). Our calculations for AGN jets show 
similar levels. A detailed discussion of time-depende nt circu- 
lar po larisation in microquasars is presented in 1 Saxton et al.l 
(|2010l l and a time-depedent polarimetry of AGN jets will be 
presented elsewhere (Saxton et al., in preparation). 



4 CONCLUSIONS 

We performed hydrodynamic simulations for jets encoun- 
tering external media with conditions and density struc- 
tures appropriate for active galaxies, and carried out time- 
dependent radiative transfer calculations to determine their 



polarization and emission properties. Our polarised radia- 
tive transfer formulation takes account of emission, absorp- 
tion, re-emission, Faraday rotation and Faraday conversion. 
The radiative transfer equations were solved explicitly fol- 
lowing the jet evolution. We applied this method to model 
the temporal evolution of emission from AGN jets interact- 
ing with a variety of structured ambient media. Our cal- 
culations showed that jet emissions vary considerably even 
though the launching conditions remain steady at the nu- 
cleus. Their variations are affected by the external density 
profile; thus environmental factors play an important role in 
determining the observed morphology and temporal proper- 
ties of radio jets. 

Our simulations show that the ambient media influence 
the distribution of jet knots and the relation between emis- 
sion intensity and advance of the jet. Determining the jet 
knot distributions and the relation between intensity emis- 
sion intensity and growth of jets can constrain the proper- 
ties of the background gas. For instance, in an effectively 
closed cocoon (that contributes much of the radio lumi- 
nosity) asymptotic power-laws appear to relate the elapsed 
time, jet length and total intensity. The indices of these rela- 
tions depend on the density gradient of the ambient medium. 
For models with an open left boundary (with the jet termi- 
nus far from the nucleus) the radio intensity appears to flash 
independently of the jet length. The total intensity can fluc- 
tuate by factors of a few, and the hotspot can fluctuate by 
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Figure 11. Same as Figure [Jl but for the fastIrc jet simulations. In the simulations, a fast jet propagates in a uniform background 
medium, and a closed boundary condition is adopted. The jet plasma is relativistic, with 7 = 4/3. 
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Figure 12. Synchrotron radio lightcurves (at 6 cm) of the fast2o jets. In the simulations, the jet traverses a background medium with 
a declining density. The left boundary is open. The magnetic field is random or ordered {B \\ v) with (3 values as indicated in each 
panel. The open left boundary reduces accumulation of the backfiow, and he jet shocks outshine the cocoon emission, (cf. the fastIo jet 
simulations). The fiuctuations of the emission from the hotspot are visible in the /hs lightcurves. 
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Figure 13. Synchrotron radio lightcurves of the FAST2 jets, as in Figure fT2] The simulations adopt same conditions as those in the 
fast2o jet simulations, except that the left boundary is refiective. 




Figure 14. Synchrotron radio lightcurves (at 6 cm) of the FAST3 jets. In the simulations, a fast jet propagates in a medium with a 
rising density. The jet generates an opaque fireball with no identifiable hotspot. 
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Figure 16. Evolution of intensity along the midline (projection of the jet axis) for the fast jet simulations (as annotated in respective 
panels, and assuming (3 = 10). In cases with an open boundary condition (fastIo, fast2o) the jet shocks are brighter relative to the 
cocoon, and the hotspot is brighter relative to the pinch shocks. The density ripples in the background medium make little difference: fast4 
resembles fastI. A radially declining density profile enables the jet's advance to accelerate (fast2, fast2o) while a radially increasing 
profile restrains the expansion of a decelerating fireball (fast3). Except in FAST3, the first jet shock is brightly visible, fluttering about 
its mean position. These flgures are derived from Stokes maps made with ordered flelds, and each vertical stripe is separately normalised. 



tens, within dynamically brief intervals (though longer than 
a human lifetime). Caution is needed when inferring the jet 
power from a single-epoch observations. In any particular 
epoch, the two hotspot s may differ by three or more times 
in intrinsic brightness, even given identical and constant jet 
fluxes, and neglecting orientation and beaming effects. 
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Figure 18. Illustration of the variability of the luminosity ratio of two hotspots in open-boundary simulations. (fastIo, fast2o and 
fastIro from left to right respectively). The magnetic fields are quasi-poloidal, B \\ v. The top row shows geometric mean hotspot 
ratios obtained from the structure function of ln(/hs), for frames when the bow shock has advanced to Xg ^ 60 kpc. The bottom row 
shows cumulative distribution functions of the hotspot ratios at independent times. The dotted, grey and dashed curves depict cases 
with /3 = 1, 10 and 100. 
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APPENDIX A: AN REMARK ON RADIATIVE 
TRANSFER 

In our calculations, radiative transfer is computed in eule- 
rian cells that do not co-move with jet plasma. Thus, the 
transfer calculation has not included relative Doppler shift 
explicitly. We now assess whether or not the omission af- 
fects the calculated results significantly. The quantity 
is Lorentz invariant (where ly is the frequency and lu is the 
specific intensity). The radiative transfer equation in a co- 
moving frame thus takes the form: 



K 



^] m (Al) 



(see e.g. lMihalaslll970l : |Peraiahll2002l ). Here, variables evalu- 
ated at the local rest frame are denoted by the subscript 0. 
The co-moving transfer equation implies that Doppler effect s 
are of the first order in [1 — (n • v)/c] (see also ICast"Qrlll972l l , 
where v is the medium local velocity. If ignoring terms of 
0[{v/c)'^] or higher and performing a series expansion, then 
the inclusion of Doppler shift simply modifies the propaga- 
tion operator: 



where 



,«3{(n.V)(Ii^)} 



(A2) 



(A3) 



The [1 — (n • v)/c] term in front of the operator ds re-scales 
time and length along the ray and is unrelated to the ab- 
sorption and re-emission of the radiation, and the Faraday 
effects on the Stokes parameters. The term Q, which is pro- 
portional to Azy/z/Q (= (u — uo) / uo) , induces a "diffusive 
drift" of power across the frequency space. It can be omit- 
ted, provided that there is no sharp relativist ic velocity gra- 
dient across the computational cells (e.g. in the presence of 
relativistic turbulence). Thus, omission of Doppler effects 
might distort the perceived time scale and length of the sys- 
tems but would not affect the qualitative results obtained 
for the global morphology and polarization level of the jets. 
Simulations with explicit inclusion of Doppler effects require 
substantial modifications in the numerical algorithm and ad- 
ditional constraints set by the relativistic physics. We will 
leave this for a separate study. 



